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Abstract 

The mass of the bottom quark is analyzed in the context of QCD 
finite energy sum rules. In contrast to the conventional approach, 
we use a large momentum expansion of the QCD correlator includ- 
ing terms to order a1{m1/q^)^ with the upsilon resonances from 
e^e~ annihilation data as main input. A stable result mf){mii) — 
(4.19 ± 0.05) GeV for the bottom quark mass is obtained. This re- 
sult agrees with the independent calculations based on the inverse 
moment analysis. 
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1 Introduction 



With the Standard Model being so well established the major theoretical ef- 
fort of today's theoretical and experimental investigations goes into a precise 
determination of the parameters of the model. An important example is the 
determination of the CKM matrix elements and their relative phases in current 
B-physics experiments. A better knowledge of these parameters will lead to a 
better understanding of the structure of weak interaction including CP viola- 
tion, and possibly the discovery of " new physics" . For an analysis of most of 
the experiments in the B-sector a precise value of the bottom quark mass is 
essential. 

In the last few years, there has been increased activity with the aim of de- 
termining the bottom quark mass [1 — 10] to higher accuracy. The most popular 
method was the one based global duality in QCD. In so-called inverse moment 
sum rules derivatives of the two point correlation function of the b quark vector 
current are compared with experimental information from electron-positron an- 
nihilation into hadrons 121 ■ These sum rules seem to be a suitable method since 
it is very sensitive to the heavy (b-quark) quark mass. Recent analysis using 
inverse moments involving 10 to 20 inverse powers of the squared energy can be 
found in |2,^ and, using a lower power inverse moments, in |^. Non relativistic 
QCD, potential models (for a review see ^31) and lattice techniques have 
also been considered in order to determine the heavy quark masses. However, 
in the case of the lattice techniques and potential models, there exist problems 
with the choice of the mass definition (running-, pole-, threshold mass) and the 
size of non-perturbative effects. So far, the available experimental information 
relies mainly on the knowledge of the masses and widths of six upsilon reso- 
nances, three of them lying below the continuum threshold of the BB meson 
production. 

The method we propose here is somewhat orthogonal and complementary to 
the one based on the inverse moment sum rules: we employ on the theoretical 
side the large momentum expansion of QCD, i.e. an expansion in m^/s, where 
s is the square of the CM energy. Such an expansion makes sense as long 
as s is far enough away from the continuum threshold. We will consider the 
perturbative expansion up to second order in the strong coupling constant and 
twelve powers of the b-mass over energy ratio using the results of the reference 
|14| . An expansion to the third order in the strong coupling but with only four 
powers in the mass-energy ratio is also known [15 | , which we will not make use 
of for reasons of consistency. On the phenomenological side of the sum rule 
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we will consider the six upsilon resonances. In addition, our method should be 
sensitive to the poorly known continuum data. We circumvent this problem 
by a judicious use of quark-hadron duality. Above the resonance region, where 
experimental data are poorly known, we incorporate a real polynomial in the 
sum rule or, what is the same, a suitable combination of positive moments, in 
such a way that the contribution of the data in the region above the resonances 
can be practically eliminated. The method has been successfully checked in 
the charm sector ^7], where the charm quark mass was predicted by using the 
experimental data and the polynomial insertion in the intermediate region. The 
result for the charm mass was found in good agreement with the ones obtained 
using other methods. Of course, in the b-sector one is in a truly heavy quark 
regime and it is not clear that the same approach may lead to a result for the 
b-mass of competitive accuracy. We will show that this is actually the case. 

The plan of this note is the following: in the second section we briefly review 
the theoretical method used here which is based on weighted finite energy sum 
rules, in the third section we discuss the theoretical and experimental inputs 
used in the calculation and present our results for the bottom quark mass with 
a discussion of the errors. We finish the paper giving the conclusions. 

2 The calculation 

In order to write down the sum rule relevant to our case, we apply Cauchy's 
theorem to the two point correlation function n(s), with the appropriate flavor 
content in the quark currents, and weight the integration with a polynomial 
P{s). The integration path extends along a circle of radius |s| = s^, and along 
both sides of the physical cut [sq, s/j]. The polynomial P{s) does not change the 
analytical properties of n(s) in the integration region, so we obtain the following 
sum rule 

- [ P{s)lmU{s) ds = / P{s) n(s) ds, (1) 

J so 27ri J\s\=sn. 

where sq is the physical threshold of the bb channel, starting at the first bb 
resonance below the continuum threshold sc ^ sq. The integration radius sr 
is chosen in such a way that on the circle the asymptotic expansion of QCD, 
IIqcd(s), constitutes a good approximation to the two point correlator. The 
left hand side of equation is related to the experimental e'^e^ annihilation 
cross section via the unitarity relation, 

R{s; e+e- hadrons) = 12n ^ Qj Imn(s). (2) 

flavors 

In order to perform analytically the contour integration in we take, as a 
matter of convenience, the two point correlation fimction nQCD('?^) which has 
been calculated in (M) as an expansion up to second order (three loops) in the 
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strong coupling constant as and as a power series of m} jq^ up to the sixth 
order, 

nQCD(<z^) = n(«)(g2) + n(i)(g^) + (^^)'n(^)(g2), (3) 

where the different terms of the expansion in ols have the form: 

n«(<z^) = EE4iK^)(i-^)'(^)'- (4) 

j=0 fc=0 ^ / \H / 

In equations H3I4(I /i is the renormahzation scale of the perturbative calculation, 
m is the running quark mass and Aj*],(m, ^) are the coefhcients of the QCD 
perturbative expansion which depend on m and /i through powers of Inm^Z/i^ 
|14| . In the asymptotic expansion of equation Q one should include the non- 
perturbative terms due to the quark and gluon condensates, which are known 
up to order |19| . Their contribution, however, is completely negligible in the 
range of energies of interest to us. For this reason we will discard them from 
now on. 

As commented in the introduction, the experimental side of the sum rule 
is dominated by the six well established bb upsilon resonances. With this ex- 
perimental information in the narrow width approximation we have the hadron 
cross section of equation (01 given by 

^(*) = E ^rcs 5{S - M2 J, 

Oi ^ — ^ 
em j.(,g 

where Mrcs and Frcs are respectively the masses and electronic widths of the six 
resonances, and aem = 1/131.8 is the electromagnetic coupling constant taken 
at the typical scale of 10 GeV, where the resonances are produced. 

Furthermore, in the theoretical side of the sum rule, equation P(s) is 
taken to be a third degree polynomial 

P{s) = ao + flis + a2.s^ + a^s^ , 

whose coefficients are fixed by imposing a normalization condition P(so) = 1, 
and requiring that it should minimize the contribution of the continuum [sc , s p;\ 
in a least square sense, i.e., 

I s"P(s) ds = for 71 = 0,1,2, 

J sc 

where sc is the value of the continuum physical threshold for the BB meson 
production. The choice of this polynomial guarantees that the contribution of 
the experimental data in the continuum region (which is experimentally poorly 



4 



known) will be negligible in the sum rule as compared to the contribution of the 
lower resonances below the continuum threshold. In fact, a smooth continuum 
contribution which can be described by a second order polynomial vanishes 
identically. In this way, we avoid the inaccurate experimental information in 
the bb continuum region. This approximation procedure was checked in the cc 
channel, where there is good experimental information on the continuum due 
to recent results from the BES II collaboration |[18j. Accurate and consistent 
results for the charm mass have been obtained, by either incorporating the 
continuum or eliminating it by the procedure above |17j . 

Notice also that the dependence on the method, which relies on the choice 
of the polynomial is very weak. To check that fact we have performed 
the same calculation using polynomials of different degrees as well as changing 
the boundary conditions of the third degree polynomial although keeping 
the same philosophy as commented in the paragraph above. For instance for a 
two degree polynomial we find a difference of .01 Gev in the bottom mass but 
the experimental error coming from the uncertainties in the resonance data is 
double than the one we find with a degree three polynomial. On the other hand 
a polynomial of three degree changing slightly the boundary conditions gives 
just the same result. 

Higher order polynomials, as well as the inclusion of the incomplete order 
in equation (PJ can be considered as well. However, before doing so, a detailed 
study of the different moments involved in equation should be performed. 
We defer this study to a further publication. 

With our study we conclude that our selection for the polynomial do not 
preclude the results found nor the errors presented in this work. 

The integrals that we have to calculate on the right hand side of the sum 
rule, equation J^l, are 

J{p,k) = ^(f sp(\n^] ds, 

for A; = 0, 1, 2, 3 and p = —6, —5, .., 3. Their analytic evaluation can be found in 
reference [201. After integration, equation ^ can be cast into the form 

97r(131.8)2^M„. r,e. PiMl^) 

res 

ri=0 i=0 j=0 k=0 

giving a non-linear equation with the quark mass (m) as the only unknown. To 
proceed further and solve this equation we have to choose a suitable value of 
sn for which perturbative QCD gives a good approximation for the correlation 
function. Furthermore, the scale /i can be chosen arbitrarily to give a prediction 
of the quark running mass at this scale. For calculational convenience, we take 
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/i^ — sr since most of the logarithmic terms vanish after integration. Then, 
after solving equation (jSJ, we run the quark mass from that scale to the mass 
scale itself, m(yLt = m), by means of the corresponding renormalization group 
equations at the appropriate loop order in the strong coupling constant j21ll22| . 

As far as the theoretical input matches the experimental data within the 
error bars, the results should be independent of the choice for sr. Nevertheless, 
the lack of experimental data in the continuum region of the bb channel makes 
sensible in our method to study the sr dependence of the quark mass results. 
Although this dependence is going to be very small, due to the choice of our 
polynomial in the sum rule, we will take the quark mass prediction at the most 
stable value with sr. 



3 Results 

The experimental inputs in our calculation are as follows: 

Firstly, the physical threshold sq is the squared mass of the first low lying 
resonances in the bb channel, T(9460), 



So = rriY 



= 9.46^ GeV^ 



whereas the continuum threshold sc is taken at the squared energy of the BB 
meson production 

sc = (2mB)2 = (10.58)2 GeV^ 

Secondly, the absorptive part of the two point correlation function is ob- 
tained from the six known bb vector resonances T(1S'),..,T(65') with the follow- 
ing masses and electromagnetic widths j23| |24|. 

M{1S) = (9460.30 ± 0.26) MeV, r(lS') = (1.32 ± 0.05) keV, 

M{2S) = (10023.26 ± 0.31) MeV, r(2S') = (0.520 ± 0.032) keV, 

M{3S) = (10355.2 ± 0.5) MeV, r(3S') = (0.48 ± 0.04) keV, 

M{4:S) = (10580.0 ± 3.5) MeV, r(4S') = (0.248 ± 0.031) keV, 

M{5S) = (10865 ± 8) MeV, r(5S') = (0.31 ± 0.07) keV, 

M{6S) = (11019 ± 8) MeV, r(6S') = (0.130 ± 0.030) keV. 



(6) 



Finally, for the strong coupling constant we take the input value at the mass 
of the electroweak Z boson |23j 

a^(A^z) = 0.118 ±0.003 (7) 

and run down to the scale of the contour radius sr, using the four loop formulae 
of reference ■ 

Now we proceed as follows. For fi^ — sr with sr, in the range [150, 350] GeV^ 
we determine m^^ (/i) by solving equation ^ keeping terms in the perturbative 
expansion up to order (os/tt)* for i = 0,1,2. Then we run the results from 



6 



= sr to fi = rUb with the appropriate renormahzation group equations 
Finally, we choose the value of sr which is most stable in the range of energies 
considered. In this way we find stability at sr = 240 GeV^ and the following 
results for the bottom mass: 

4.09 GeV, 

4.23 GeV, (8) 
4.19 GeV. 

To estimate the various errors arising in the calculation of the b quark mass 
we consider first the uncertainties in the masses and widths of the resonances in 
equation ©. This gives an experimental error for the mass Eexp = 0.02 GeV^. 
Secondly, we consider the uncertainty in the strong coupling constant, equation 
Q, which leads to Sa = 0.01 GeV. Finally, we include a conservative asymp- 
totic error originating from the higher orders in expansion of equation © . To 
estimate this error we consider the difference between the second and first order 
results which gives easy = 0.04 GeV. Adding the errors quadratically, we find 
for the b-quark mass 

nibimb) = (4.19 ± 0.05) GeV. (9) 

Before going on, a comment on the non perturbative contributions to the two 
point correlation function is in order. Using accepted values for the condensates, 
we find that their contribution to the b-quark mass for the considered range of 
sr is always of the order of 0.001 GeV and, therefore, negligible in comparison 
to the errors obtained before. We have also checked that the convergence of the 
result © with respect to the power series expansion in m^/s has no influence 
to the error bars in the sr stability domain. 

To have a flavor of the whole procedure, we include two figures which sum- 
marize the main steps followed in the calculation of the bottom quark mass. 

As a sample calculation, we have plotted in figurenthe results of ml^'^\mb) 
as function of the radius sr. The upper (lower) curve corresponds to the first 
(second) order in the strong couphng constant. We find the most stable result 
for Sr = 240 GeV^ . The difference between second and first order results is 
Easy — 0.04 GeV. The sr dependence is so tiny that all the values for the mass 
in the whole range [150, 350] GeV^ lie within the relatively small asymptotic 
error bar. 

In figure El we have plotted the right hand side of the mass equation JSJ as a 
function of mb{fJ.) with the renormahzation point taken at the value of sr and 
the latter in the stability point (/Lt^ = sr = 240 GeV^). The thin curve includes 
only the calculation to one loop, whereas the medium and thick curves include 
calculations to two and three loops, respectively. The straight strip corresponds 
to the contribution of the experimental data to the left hand side of equation ((SJ 

■^Recently, it has been claimed an enhancement of the experimental data in the 45 and 5S 
resonance region 1251. This could change the final value of the bottom mass by less than 0.01 
GeV, which is within the cxprimcntal errors of the method. 



mrimb) = 
ml^\mb) = 
m^'\mb) = 
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Figure 1: Stability of the two (upper curve) and three (lower curve) loop calcu- 
lations of the bottom quark mass. 

including experimental error bars. Hence, the crossing of the straight strip with 
the three curves gives the solution of equation ||SJ), i.e. the results for the running 
mass m^°'^'^^(240 GeV^). Running this results to fi — mi, we find the values 
quoted in equation |(SJ). Notice the very good convergence of the perturbative 
series at the energies relevant to this process. 

4 Conclusions 

The bottom quarks mass has been determined in an unconventional fashion by 
employing finite energy QCD sum rules with positive moments. By apply- 
ing Cauchy's theorem to the bb vector current correlator weighted with a real 
polynomial, it was possible to virtually eliminate the contribution of the yet 
unknown experimental data in the continuum region from just above the res- 
onances to the beginning of the asymptotic region where QCD is valid. The 
experimental input needed is only the resonance masses and couplings. 

The method we use is quite independent of the more popular one based on 
inverse moment sum rules. The result we obtain, mt{mi,) — (4.19 ± 0.05) GeV, 
agrees with most of the recent calculations especially with references |S] 
and that study high and low inverse moment sum rules, respectively. The 
errors of our calculation are comparable to the ones quoted in these references. 
Such independent determinations of the same quantity serves as a probe of the 
validity of the concept of duality which are at the basis of most calculations 
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m^^CGev) 

Figure 2: Results for one (thin), two (medium) and three (thick) loops running 
bottom masses at the renormalization and stability scale = sr = 240 GeV^. 

in the b-sector. We find nice convergence of our result with respect to QCD 
perturbation theory and negligible non-perturbative contributions. 
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